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ABSTRACT 

V 1  f 

'  An  algorithm  is  proposed  for  solving  one^ dimensional  free  boundary 
problems  with  change  of  phase.  The  technique  consists  of  solving  the  heat 
equation  in  progressively  increasing  rectangles  whose  size  is  controlled  by 
the  Stefan  condition.  Convergence  of  the  scheme  is  shown  and  an  estimate  of 
the  rate  of  convergence  is  given  • 
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SIGNIFICANCE  AND  EXPLANATION 


Stefan  type  problems  arise  as  descriptions  of  phenomena  such  as  melting 
of  metals,  solidification  of  alloys,  crystal  growth,  permafrost  behavior.  Of 
particular  importance  in  such  problems  is  the  shape  and  evolution  of  the  free 
boundary  or  interface  • 

This  paper  suggests  a  method  of  constructing  the  free  boundary  by  solving 
the  heat  equation  in  a  sequence  of  increasing  rectangles  •  The  interface  is 
then  approximated  by  piecewise  vertical  segments. 

The  simple  geometry,  and  boundary  conditions  suggested,  can  be  used  to 
perform  efficient  numerical  calculations. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


AN  ALGORITHM  FOR  THE  ONE-PHASE  STEFAN  PROBLEM 
E •  DiBenedetto*1^  and  R.  Spigler*2^ 

1.  Introduction; 

Consider  the  following  one-phase  one- dimensional  Stefan  problem 
La  =  uxx  -  *  0  in  =  {o  <  x  <  s(t)}  *  (0,T] 

u(x,0)  «  h(x),  0  <  x  <  s (0 )  «  b,  (b  >  0) 

(SP)  <  uy(0,t)  -  g(u(0,t) ,t) ,  0  <  t  <  T 

u(s(t)  ,t)  =*  0,  0  <  t  <  T 

Ux(s(t),t)  -  -s(t),  0  <  t  <  T 

where  x  ♦  h(x),  (£,t)  ♦  g(£,t)  are  given  functions  on  (0,b]  and  R  *  (0,T] 
respectively • 

Under  suitable  assumptions  on  h(*)  and  g,  (SP)  admits  a  unique  classical 
solution.  For  such  results  we  refer  to  the  survey  article  [13]  and  other  papers  given 
in  the  extensive  bibliography. 

The  aim  of  this  paper  is  to  propose  an  algorithm  to  construct  the  solution,  which 
consists  in  solving  the  heat  equation  in  progressively  increasing  rectangles,  whose 
size  is  controlled  by  the  Stefan  condition  ux(s(t),t)  *  -s(t)  . 

Such  an  algorithm  arises  as  a  natural  modification  of  Huber's  method  [15,10,1] 
and  can  be  described  in  a  simple  fashion  as  follows. 

First  the  interval  [0,T]  is  divided  in  n  intervals  of  length  0  *  T/n,  then 
for  t  e  [0,9]  we  set  Sg(t)  *  b  and  solve  the  problem 
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(Pi) 


uxx  •  ut  *■  0  in  R1  5  (o  <  *  <  b>  *  (0,6) 

u£(0,t)  -  g(u1(0,t) ,t) ,  0  <  t  <  6 

(x ,0 )  «  h(x) ,  0  <  x  <  b 

u^bjt)  -  0,  0  <  t  <  6  . 

We  compute  the  number  u\sq(9)  -  b,8)  and  determine  the  rectangle 
R2={0<x<x2»b-  u*(b,9)8>  *  (6,29)  , 

setting  Sg(t)  -  x^  for  t  e  (6,26)*  in  R2  we  solve  a  problem  similar  to  (Pi),  and 
proceed  in  this  fashion  • 

The  convergence  of  schemes  where  at  each  time  step  the  free  boundary  is 
approximated  by  a  vertical  segment  was  conjectured  by  Datzeff  [5,6] . 

A  proof  of  convergence  has  been  given  by  Fasano-Primicerio-Fontanella  [11] • 

Their  scheme  however  is  somewhat  more  complicated  than  the  one  we  propose  here,  both 
in  the  construction  of  the  sequence  of  rectangles  and  in  the  boundary  conditions  on 
*j,(j  ~  1)8  <  t  <  j9  which  are  not  homogeneous,  being  given  as  a  relationship  linking 
the  distance  t  -  (j  -  1)8  with  the  values  -  1)8)* 

Thus  the  scheme  we  have  described  has  a  two-fold  simplicity t  the  rectangular 
geometry  and  the  homogeneity  of  the  boundary  data  on  the  approximating  free  boundary. 
We  treat  the  problem  for  boundary  data  on  x  -  0  of  variational  type  since  such 
a  oonditlon  is  the  "natural"  one,  as  pointed  out  and  discussed  in  [7] .  The  method 
could  handle  the  Dirichlet  boundary  data  as  well* 

We  give  an  estimate  of  the  speed  of  convergence  which  turns  out  to  be  of  the 
order  of 

As  a  related  work  we  mention  the  methods  of  [14]  based  on  enthalpy 
considerations,  which  yield  a  slightly  poorer  rate  of  convergence,  of  the  order  of 
(8  In  0"1]  h  . 


The  methods  of  proof  are  simple  in  that  we  exploit  both  the  rectangular  geometry 
and  the  homogeneity  of  the  data  to  represent  the  approximating  solutions  by  means  of 
elementary  heat  potentials. 


Section  2  contains  the  precise  description  of  the  algorithm,  assumptions  and 
statement  of  results*  In  Sections  3  and  4  *e  produce  basic  estimates  and  prove  the 
convergence  of  the  approximating  solution,  to  the  solution  of  (SP)* 

The  error  estimate  is  given  in  Section  5*  We  conclude  the  paper  by  discussing 
some  variants  of  the  scheme* 


-4- 


numbers  u^(x^,j9)  are  non-positive  so  that  x^+1  >  x^  and  Sg(*)  is 
non- decreasing . 

On  the  domain 


We  can  now  state  our  main  result. 

Theorem  As  9  ♦  0,  Ug(x,t)  ♦  u(x,t)  uniformly  in  S  and  Sg(t)  ♦  s(t)  uniformly  in 
[0,T]  •  Moreover  there  exists  a  constant  C  depending  upon  H,  b,  G ^ ,  G2,  T  such  that 

»ufl  -  ul  <  c/5 

Is.  -  Si  <  c/5  . 

9  T 

Remark t  (i)  In  view  of  the  stability  of  (SP)  (see  [2,3])  the  Lipachitz  condition  in 
[A ^  ]  can  be  replaced  by 

[A y ] 1  x  ♦  h(x)  is  essentially  bounded  in  [0,b] * 

(ii)  The  signum  condition  on  g  in  [A2J  can  be  dropped,  provided  we  assume 
g<0,t)  *  0,  t  9  (0,T]  •  In  this  case  we  set  G2  *  0  In  the  growth  condition  for  g* 

-5- 


Mu  • 


3.  Some  basic  estimates: 


4 

i\ 


r.  i»_- ji!i 

L  4(t  -  T)J 


r<x#t;?#T)  »  -  exp 

2/w(t  -  T) 


be  the  fundamental  eolation  of  the  heat  equation  and  let  G(x,t;€,T),  N{x,t;£,x)  be 
the  Green's  and  Neumann's  functions  respectively,  defined  by 


G(x,t|€,X)  a  ^(X|t|C|T)  “  ^*(“X,tj£,T)  f 
N(x,t;£,T)  »  r(x,t|£,T)  +  r(-x,t»C,T)  . 


In  the  rectangle  Rj  the  solution  u^  of  (P^)  can  be  implicitly  represented  as 

xj 

(3.D  u3(x,t>  -  /  N(x,t>e,(j  -  i)9)uj<e,(j  -  n@)d«  - 


/  N(X,t»0,T)g(u3(0,T) ,T)dT  + 

( j-1 )0 

t  . 

/  N(x,t;x. ,t)u3(x  ,T)dT 


(j-1)0 


j  -x-  j- 


for  (x,t)  e  R j  •  Taking  the  derivative  with  respect  to  x  in  (3*1)  and  letting 


x  ♦  x^  we  obtain 


Xj 

(3*2)  ~u^(x.,t)  *  /  G(x  ,t;S,(j  -  1)9)u^(5,(j  -  1 )  © ) )  < 

«  x  J  n  1  * 


-  /  H  (X  ,t|0,T)g(u3(0,T) ,T)dT  + 

(j-1 >o  3 


t  . 

*  /  N  (X.,t>X,,T)ui(x,,T)dT,  (j  -  1)0  <  t  <  j0  . 

(1-1  )0  *  3  3  X  3 


The  calculations  leading  to  (3#1)-(3*2)  are  routine  and  we  refer  to  [4,9,12]  for 


details  • 
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I 


I 


Let  us  fix  1  <  j  <  n 


over  Rk,  1  <  k  <  j.  Since 


and  (x,t) 

h  (Nuc  ■ 

in  Rk  we 


ft  Rj  and  integrate  the  Green  identity 

d 

uN^)  -  (Nu)  *  0 

are  away  from  the  singularity  we  obtain 


(3.3) 


k  * 

N(x,t;e,k0)uK(C,k0)dS  -  /  N(x,t>£,(k  -  1)0)uk(£,(k  -  1)0)d£ 


k0 

*  -  /  N(x,t|0/T)g(Ufl(0#T)  ,T)dT 

(k-1 )6  9 


k0 


/ 

(k-1 ) 0 


N(x,t;xk#T)ux(xk/T) dT 


for  k  «  1,2, ...,(j  -  1) . 

By  virtue  of  our  definition  of  (Pj),  the  second  integral  on  the  left  hand  sid*  >f 
(3.3)  can  be  rewritten  as 


Vi 

-/ 


N(x,tiC,(k  -  DSju^1  (£,(k  -  1 )  0 )  d£  . 


Consequently  adding  the  identities  (3.3)  for  k  -  1,2,...,(j  -  1)  with  (3.1)  and 
recalling  the  definition  of  uQ(x,t)  we  obtain 


b  t 

UQ(x,t)  *  /  N(x,tj£,0)h(£)d£  -  /  N(X,tf0,T)g(Uft(0,T)  ,T)dT 

0  0  9 

j-1  k0 

+  I  /  N(x,t;x  ,T)uk(x  ,T)dT  + 

k-1  (k-1 ) 0  x  x  K 

t 

+  /  N(x,t;x  ,T)u^(x.,T)dT  . 

(j-1 ) 0  *  x  * 


Lemma  3.1:  For  each  0  the  following  estimates  are  valid 


(a) 


2G  /T  2G-ST/S* 

lu0<°^>lT  <  fHb  +  1®  =  CQ  , 


( b> 


< 


*,s 
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Proof i  By  virtue  of  the  maximum  principle  ufl  >  0  in  V  and  uj(x4,t)  <  0, 

a  se  3 

t  e  ((j  -  therefore  dropping  the  non-positive  terms  on  the  right  hand  side 

of  (3.4)  and  letting  x  +  0  we  obtain 


b  .  t  g(u0(O,T),T) 

0  <  ufl(0,t)  <  /  2r<0,ti5,0)h(£)d£  -  —  /  - -  ■■ - 

"0  /w  0 


dt 


b  .  0  <5  t  ufl(0,T) 

<  —  /  —  expl-£2/4t]  dC  +  —  /  ■■ 

/*  0  /t  /wo  /t  -  T 


/t  -  T 
2G_ 

dT  +  -  /T  . 

/w 


Statement  (a)  is  now  a  oonseguenoe  of  standard  calculations  and  Gronwall's 
inequality • 

Statement  (b)  follows  from  the  maximum  principle  applied  recursively  to  (P^). 
Lemma  3,2:  For  each  j  *  1,2,...,n 

(a)  |u£(Xj,t)|  <  HeT  +  8G*e*  iC,,  *  -  T//^b2 

3 

(b)  l*J^  Uq (x,t)  |  <c1#  (x,t)  e  Pa  , 

s 

where  G  -  GjCq  +  G2,  and  H  -  max{H,G> . 

Proof?  We  employ  an  induction  argument  by  making  use  of  formulae  (3.2),  First  we 
prove  that  if  for  some  j  -  1,2, ...,(n  -  1)  we  have 

I*£(x,j0)  |  <  P,  0  <  x  <  Xj 
for  some  positive  constant  P»  then 

(3.5)  l'4+1(xi+1't)  I  4  Pe9^15  +  8G  -2—  e9/,^b  ,  t  e  ( j0,(  j  +  1)8) 

/»b 


where 


max  |g(u0(O,t)  ,t)  |  <  G^  +  G2  s  G  • 
( 0  ,T]  * 


Consider  (3.2)  written  for  the  Integer  j  +  1 

(3.6)  |  uj+1(xj+1,t)  -  X,  +  I2  +  I3 

and  estimate  the  integrals  1^,  i  *  1,2,3  separately  as  follows. 


i 


-8- 


j+1 


I1I  *  l/  G(x 4+1 't;€,  j9)u^  ( C ^  j 0 >  aC I  *  by  definition  of  (P^  ) 

0  3 

Xj  Xj 

-  j/  G(x  ,tjC#j6)u^(C,j6)dC|  <  P  /  lG(x  #t;5,j9)|dC 


j+V 


* 


j  +  1‘ 


*  P  r  -n  P 

<  —  /  e  an  -  - 

/wo 


To  estimate  13#  I3  we  recall  the  following  elementary  estimates  on  Nx 


'Wi't,0'T)'  4 


/wb‘ 
1 


2  ' 


(xj+,  >  b  >  0) 


Therefore 


,I2I  «  (t  -  j0)  <  40  T~i  ' 

fnb  /*b 


Putting  together  these  estimates  as  parts  of  (3.6)  we  have 


for  all  t  e  (j9,(j  ♦  1)6). 

Consequently  by  Gronwall's  inequality  (3.5)  follows  at  once. 

Consider  now  the  problem  ( P  ^ )  j  ■  1,  Since  x  ♦  h(x)  is  Lipschitz  continuous 

in  [0,b]  ,  h*(x)  exists  for  a  .e .  x  e  [0,b]  and 

ess  sup  |h*(x)|  <  H  <  H  -  max  (h,g)  . 

t0,b] 

Therefore  by  the  previous  argument 
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|u’(x1ft)|  <  He*  +  0G3®a#  t  e  (0,6] 
where  for  simplicity  of  notation  we  have  set  a  •  8//wb2  .  The  function 
(x,t)  ♦  u^(x,t)  5  v(x,t),  will  satisfy  the  Dirichlet  problem 
Lv  -  0  in  R1 

v(0,t)  -  g(Ug(0,t> ,t)/  |g(uQ(0,t)  ,t)  |  <  G,  te  (0,6) 

|v(x1#t)|  <  He*  +  8G<*e*,  t  e  (0,0] 

v(x,0 )  -  h'(x),  l h *  ( x)  |  <  H  a.e*  x  6  (0,bj  . 

Consequently  by  the  maximum  principle 

|u\x,8)|  <  He*  +  0Gae* 


and  by  (3*5) 


|u2(x2,t)|  <  He20  ♦  8Gae2*  ♦  8Gaea 

<  H*20  ♦  8G(2<*).2<\  t(9,20]  . 


Proceeding  in  this  fashion  we  obtain 

|u£(xj,t)|  <  He^a  +  BGtjoOe^0  . 

Now  ja  -  -  ■  j9  <  T,  and  therefore  the  lemma  is  proved*  Next  we  introduce  the 

/*b2 

function  t  ♦  Bg(t)  defined  by 

(3.7)  ;  (t)  -  x.  -  u^(x.,j8)(t  -  (j  -  1)6),  te  U)  -  1)efj6)  . 

e  j  x  j 

For  t  *  jO,  ig(t)  -  xj+v  *q(0)  80  the  graph  of  a g ( • )  is  obtained  by 

connecting  the  points  (b,0 ) , (x* ,0) , . « « , (x  ,(n  -  1)0)  for  t  e  (0,(n  -  1)6],  and  by 

2  n 

the  graph  of  (3.7)  for  te  [  ( n  -  1)6,t]  .  The  points  (x^,(j  -  1)8),  j  -  1,2,...,n, 
are  the  lower  vertices  at  the  right  side  of  the  R^'s. 

By  Lemma  3*2 

b  <  Sq (t)  <  b  +  CtT,  t  e  (0 ,T] 

therefore  the  sequence  (sgt*))  is  equibounded  and  equi lipschitz ,  so  that  by  Asooli- 
Arzela  theorem  a  subsequence  relabeled  with  9,  converges  uniformly  to  some  non¬ 
decreasing,  Lipschitz  continuous  curve  t  ♦  s*(t),  with  Lipschitz  constant  bounded 
by  Cj  •  Since  lig  -  Sg>T  *  C1^f  also  Sg(t)  converges  uniformly  to  s*(t)  • 
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Let  Ds*  be  the  domain  defined  by 

Dfl#  =  {O  <  x  <  s*(t)}  *  (0,T]  , 
and  let  u*  be  the  unique  solution  of 

Lu#  *  0  in  Dg#  , 

u*<0,t)  -  g(u*(0 ,t) ,t) ,  te  (0  ,T)  , 

(P)  1 

u*(x,0)  -  h(x) ,  x  e  (0,b]  , 

u*(s(t),t)  *  0,  t  e  [0,T]  • 

we  will  show  that  ufl(x,t)  ♦  u*(x,t)  uniformly  in  S  and  that  the  pair  (u*,s*)  so 
obtained  is  actually  the  unique  solution  of  (SP)  in  the  introduction. 

We  remark  that  as  a  consequence,  in  view  of  the  uniqueness  for  (SP),  the 
selection  of  subsequences  is  superfluous. 

The  following  lemma  will  be  needed. 

Lemma  3.3;  Let  (x,t)  ♦  v(x,t)  be  the  unique  solution  of 

Lv  *  0  in  Dg*  , 

vx(0,t)  -  g(uQ(0,t) ,t) ,  te  (0,T)  , 

v(x,0)  *  h(x),  x  6  (0,b]  , 

v(s*(t) ,t)  *  0,  t  6  (0 ,T)  . 

Then 

0  <  v ( x , t )  <  C1(s*(t)  -  x),  (x,t)  e  D8#  • 

Proof ;  The  lemma  is  proved  by  standard  barrier  techniques  and  the  maximum  principle 
[101  • 
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4,  Convergence  of  the  ichfn 

Lemma  4  *  1 1  Ug(x,t)  ♦  u*(x,t)  uniformly  in  8,  as  9  ♦  0, 

Proof x  By  the  triangle  inequality 

•uA  *  u*l  <  lw. I  +  lw  1 

9  u  «,8  1  -»S  2  -,3 

where  Wj(x,t)  ■  Ug(x,t)  -  v(x,t)  and  w2(x,t)  »  v(x,t)  -  u#(x,t),  and  v  is 
defined  in  Leama  3.3. 

Set 

a(t)  ■  «in{s0(t) #s#(t) 9(t)  »  aax{»0 (t) ,s* (t)}  , 

5<t)  -  0<t)  -  a(t) ,  t  e  [0,T]  . 

We  already  know  that  5(t)  ♦  0  uniformly  in  [0,T) ,  as  B  ♦  0. 

Consider  the  rectangle  Rj,  1  *  j  <  n.  We  claim  that  if 

-  D0)l  <  x  e  then 

|wt(x,t)l  <  Ct  ■«  «t ,  (x  ft)  6  <0,**)  X  ((j  -  1)9  ,j9]  . 

If  for  t  e  t(j  -  1)9,j9],  x^  <  s*(t)  then  w1  solves  the  problem 

Lwt  -  0  in  {o  <  x  <  x^}  *  ((j  -  1)9, j0]  , 

w1  (0,t)  -  0,  t  e  ((j  -  1)9, j9j  , 

x 

Iw^Xffj  -  1)9)1  <  C1151  ( j-!  )Q'  *  e  (0'xj1  * 

w^x^t)  -  -v(Xj,t),  te  ((j  -  1)9, j9j  , 

Hence  |w,.(x,t)|  <  max(c  l$1  max  v(x  ,t)>,  (x,t)  6  Ri.  By  Lentna  3 

(j  }  l  ( j”1 )9, j6]  j  J 

obtain 

|w,(x,t)|  <  c1 *6,t ^  <*,t)  e  to,-)  *  [(j  -  i)0,j0]  . 

If  for  t  e  [tj  -  1)0,j0],  x^  >  s*(t),  then  solves  the  problem 

Lw,  -  0  in  (0  <  x  <  s*(t)}  x  ((j  -  1)6,10]  , 
w1x(0,t)  -  0,  t  e  ((j  -  1)8, j0]  , 

|w,(x,(j  -  1)0)1  <  C1,5.(  0  <  x  <  s*I(j  -  1)0]  , 

w,(s*(t),t)  -  ue(s*(t),t),  te  ((j  -  1)0,j0]  . 

By  Lemma  3 .2  we  have  0  <  u0(x,t)  <  C,  <Xj  -  x),  (x,t)  e  R^  and  therefore 

I w, ( x , t ) |  <  c,,5,t.  ( x , t )  e  (0,-)  *  ((j  -  i)0, j0]  . 


3  we 


-12- 


\.-  J* 


•1 


If  t*  e  ((j  -  1)8, j0)  such  that  x^  -  s*(t*),  then  we  repeat  analogous  arguments 
in  the  domains  so  determined. 

Now  since  for  t  -  0,  (x,0)  -  0,  x  e  (0  ,*) ,  an  inductive  arg\»ent  gives 

I w j (x, t)  |  <  c1i«it,  xe  (0,-),  te  [0,T]  . 

As  for  w2#  since  it  solves  the  problem 

I*2  “  0  in  Ds*  • 

w2x(0,t)  *  g(ug(0,t) ,t)  *  g(u*(0,t),t)  , 
w2(x,0)  -  0,  x  e  (0 ,b) 

w2(s*(t),t)  -  0,  t  e  (0 ,T] 

it  can  be  dominated  [2,3]  by  the  function 


w  (x,t)  -  /  N(x,t;0,T)|g(u8<0,T),T)  -  g(u*(0,T),T)|dT 

0 


the  unique  solution  of 


“0  in  S  , 

w2x<°,t)  -  -|g<u9(0,t),t)  -  g(u*(0,t),t)|,  t  e  (0,T]  , 


w2 (x,o )  -  o,  x  e  (0,“)  . 


|w2(x,t)|  <  Gt  /  N(x,tjO,T)|u0(O,T)  —  u* (  0  *  X  ) | dT 
0 


deduce  that 


-  t  | uQ (0 , t )  -  u#(0,x)| 

lurt(x,t)  -  U*(x,t)l  <  Iw^x^t)!  +  G„  —  /  -  dT 

6  1  1  n  o  /rrr 


sup  |Ug(x,T)  -  u(X,T) J 

+  !ij  aioia. - dt  . 

i  t  j:  *  jt~~z 


fl 


And  by  Gron  wall's  inequality 


(4.1) 


2G 

— L  fi 

■up  |uB(xft)  -  u(x,t) |  <  C  151  e  *  -  C,l$l 

xe[0,*)  1  r  4  t 


for  all  t  e  [0,T]  •  The  lemma  is  proved* 

Lemma  4.2*  The  pair  (u*,s#)  coincides  with  the  unique  solution  of  (SP)  ♦ 

Proof*  The  only  thing  that  remains  to  be  proved  is  the  Stefan  condition 

u(s(t),t)  *>  -a(t),  t  e  (0#T]  •  Such  a  condition  has  been  shown  to  be  equivalent  to  the 

Integral  Identity  [2,  page  85] 


t  b  s(t) 

(4.2)  a(t)  m  b  -  /  g(u(0 ,t) ,t) dT  ♦  J  h(x)dx  -  /  u(x,t)dx  , 

0  0  0 


and  hence  it  will  be  sufficient  to  prove  that  s*,  u*#  satisfy  (4.2). 
Integrating  the  equation  Lu^  -  0  over  Rj  we  obtain 

j* 


-  /  u^(x  ,T)dT  -  -  /  g(u^  (0,T) #T)dt 


(j-t)9 


*  j 


(4.3) 


j  j 

+  J  u^(x,(j  -  1)0)dx  -  /  u^(x#jS)dx  • 
0  0 


Also  for  (x,t)®Rp*  1  <p<n#  Integrate  Lup  ■  0  over  the  rectangle 

(0  <  x  <  x  }  *  ( (p  -  1 )0,t]  .  It  gives 
P 

.  .  x 

t  t  p 

(4.4)  -  /  up(x  ,T)dT  -  -  /  g(up(0,T) ,T)dT  +  /  uP(xf(p  -  1)0)  dx 

<P-1)0  x  P  (p-1 )  0  0 


-  /  up(x,t)dx  • 
0 


By  the  definition  of  (P^)  we  have 


,  1*1  ,  1  i 

/  (x,j0)dx  -  /  u3(x,j0)dx  , 

0  0 
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therefore  adding  the  identities  (4*3)  for  j  ■  -  1  and  (4.4)  we  obtain 


p-1  j0  t 

-  I  J  uj(x  #T)dT  -  /  u^(x  ,t)dT 

j-1  (j-t)d  x  3  (p-1)9  X  P 


t  b  p 

/  g(ufl(0,T) ,t)dt  +  /  h(x)  dx  -  /  up(x,t)dx  • 
0  0  0 


We  rewrite  the  left  hand  aide  of  (4.5)  aa  f o 1 Iowa : 


p-1  j0  t  p-1  j0 

I  /  UMX ,,T)dT  -  /  U^(X  ,T)dT  -  -  l  I  <(X.,j6)dT 

j-1  ( j-1 )0  *  3  (p-1 ) 8  p  j-1  (j-1 )6  3 

t  p-1  j9  .  t 

-  /  up(x  ,p9)dx  -  l  /  e?(T) dt  -  /  ep(T)dT  , 

(p-ue  x  p  j-i  (j-i )0  (p-i)9 


e*(t)  -  uhx-.t)  -  U1(x4,i0),  (i  -  1)9  <  t  <  i8  ,  i  -  1,2, ...,n  . 

w  X  i  XI 

We  observe  that  the  numbers  -u^(x  ,j6)  are  the  slopes  of  the  Lipschitz 

x  j 

continuous  polygonal  t  ♦  s«(t)  for  (j  -  1)9  <  t  <  j9,  consequently 


p;i  ,J9 


l  J  U^(X  ,J0)dT  -  /  U*(x  ,P0)dT  -  /  -rz  ig(T)dt  -  Je(t)  -  b  . 
-1  (j-1  )0  x  3  (p— 1  )9  X  P  0  d  0 


Carrying  this  in  (4.5)  gives 


t  b  p 

:0(t)  *  b  -  /  g(u0(O,T)#t)dT  -f  /  h(x)dx  -  /  Ug(x,t)dx  + 
0  0  0 

p-i  j9  ,  t 

+  I  /  ei(T)dT  +  /  e?(T)dT  . 

J-1  (J-1)0  (p-1 ) 9 


By  virtue  of  Lemma  4.1  and  the  uniform  convergence  s0(t)  ♦  s*(t),  letting 
9  ♦  0  in  (4.6)  gives 


4 

i 

m 


V 

* 

f-.vy... 

(4.7) 


t  b  .Mt) 

s*(t)  «  b  -  J  g(u*(0 ,T) ,t) dT  +  /  h(x)dx  -  J  u*(x,t)dx  + 

0  0  0 


p-1  j0  t 

♦  li»  {  [  /  .2(T)<JT  +  /  «|(T)dT}  . 

8*o  j-i  <3-1 ) e  (p-i ) 0 


Therefore  the  lemma  will  be  proved  tf  we  show  that  the  limit  in  (4.7)  is  aero. 

In  order  to  estimate  the  «g(M  we  will  need  a  representation  for  u£(Xj,t), 
t  e  ((j  -  1)9, jO] . 

Consider  identity  (3.3).  By  taking  the  derivative  with  respect  to  x  and 


Integrating  by  parts  the  first  two  Integrals  (the  identity  Nx  *  -G^  is  used)  we 
obtain 


\  *k 

(4.8)  /  G(x,t|C,k0)uk(€,k0)dC  -  /  G(x,tf€* (k  -  1)0)uk(C,(k  -  1)0)d£ 

0  X  0  X 


k0  k0 

-  /  N  (x,tj0,T)g(ufl(0,T)  ,T)dT  +  J  N  (x,t|X^,T)uNx^,T)dT 

f  i _ 4  \  a  *  ™  _ «  t  a  "  XXX 


(k-1 )0 


(k-1 )0 


Now  by  the  construction  of  u  ,  the  second  integral  in  (4.8)  can  be  rewritten  as 


Vi 


-  /  G(x,t j€, (k  -  1)0)uk“1 (C,(k  -  1 )0) d£  , 
0  x 


therefore  adding  the  identities  in  (4.8)  for  k  -  -  1),  changing  the  sign 

and  evaluating  the  sum  at  x  -  Xj,  gives 


b  J-1 

(4.9)  J  G(x  ,t;C,0)h’ (5)<K  -  J  G(x  -  1)0)u3_1(5,<j  -  1)8)d5 

0  3  0  3  x 


-  J 

0 


( j-i )0 


Nx(x^,tjO,t)g(u0(O,T),T)dT  - 


j-1  k© 

-  I  I  N  (x  ,t;x  ,r)u  (x  ,x)dt  . 
k-1  (k-1 )6  3  *  x  * 
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Consider  now  the  representation  (3*2)  •  By  the  way  has  been  constructed,  the 

first  integral  on  the  right  hand  side  of  (3.2)  reads 


,  1-1  i-1 

/  g(x  ,t»c,(j  -  i)9)U;  a, a  -  i)0)<K  . 

0  3  x 


Finally  adding  (3.2)  and  (4.9)  gives 


,  a  b  t 

(4.10)  ?>£<*  ,t)  -  /  G(x<,t»5,0)h*  (?)dC  -  /  NJx4,ti0,T)g(uB(0,T),T)dT  + 


x'T 


j-1  k0 

+  I  /  N  (X  ,t|X ,T)UX(*  ,T)dT  + 

k-1  (k-1) 9  3  K  x  * 


+  /  H  (X  ,t»X  ,T)u^(x  ,T)dT,  te  (  ( j  -  1)9,  JO] 

(j-1 )9  3  3  x  3 


The  error  terms  e^(t)  J  “  2,3,...,(p  -  1),  can  now  be  written  as 


■je^(t)  - -j  (u^(xjft)  -  uJ(xjrj«)J  -  /  (G(x^ ,t;C,0)  -  Gfx^.jOiC.OJIh^OdC  - 

t 

"  !  (N  (X  ,t»0,T)  -  M  (X . ,J0>0,T)] g(uB (0,T),T) dT  + 

0  *  J  x  j  ° 


♦  /  M  (x.,j9»0,T)g(u  (0,T),T)dT 

t  x  j  9 

j-1  k9 

+  I  /  [H  (X  ,t»X  ,T)  -  M  (X  ,j9,X  ,T))uNx  ,T)dT 

k-1  (k-1  )9  x  3  K  X3  K  xx 

t 

♦  /  (M  (X  ,t>X  ,T)  -  M  (X  ,j9,x.,T)]u^(X.,T)dT  - 

(j-1  )9  XJJ  x  J  j  x  j 


j0 

/  Nx(xj,j0,Xj,T)u^(Xj,T)dT 


J1  +  J2  ♦  +  ^4  ♦  J4  • 
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By  the  mean  value  theorem  and  standard  estimates  on  the  exponential  function  we  have 


Ig(x  ,t»e,o)  -  g(x  ,je>e,o)  i  <  (Kl  +-372)8  <  k,9  +  k2<j  -  i)_3/2e"1/2  . 


Here  and  in  what  follows  denote  constants  independent  of  j  #  Consequently 

|J1 1  <  Hb(K^ 8  +  K2(j  -  1)“3/20“1/2)  . 

As  for  J2  analogous  arguments  yield  the  estimate 

In  (X  ,t»0,T)  -  N  (X  ,jO|0,T)|  <  K  0  , 

X  j  X  j  -J 

where  we  have  used  the  fact  that  x^  >  b  >  0.  Therefore  J2  can  be  estimated  by 

|  J2  |  <  gk3t0  * 

Let  us  now  turn  to  J3 .  First  we  note  that  for  k  -  1,2,  ••*,(;)  -  2), 

In  (x  ,tix.T)  -  n  (x  ,j0jx  ,t)|  <  (j&  -  t)  sup  |n  <x  ,pix  ,t)|  < 

X  3  X  x  3  K  xt  J  K 


a  ~  V 

<  V  *  v  ^^5  •'  410 ' 41 

4  5  t<P<J0  (p  -  t>5/2 


Now  since  (x^,k0),  (Xy(j  -  1)0)  belong  to  the  mxImI  extension  of  the  graph 
(t,8g(t)),  we  have  (Xj  -  x^)  <  JC^CP  -  t),  for  k  m  1,2,. ..,(j  -  2),  and  hence 

j-2  kO 

|J,I  <  l  I  |N  <x  ,t>x  ,t)  -  n  (x  ,10>x  ,t)||uX(jc  ,T)|dr  + 

3  k-1  <k-1)0  X  3  k  x  3  k  x  * 


<1-1)8 


+  /  |N  (X  ,t»X  ,T)  -  N  (X  ,  J0|X  ,T)  I  lu3”  (x  ,T)  I  <ax 

(j-2)0  X  j  j*1  X  J  J  1  x  3  I 


<i-2>9  <1-2)0 

<  C,K  0  /  dT  +  2C,K  0  /  [(J  -  1)0  -  t]  J/  dT  * 

1  4  0  '  9  0 


<1-1)0 

♦  c,  j  In  <X  ,t>x  ,T>  -  N  <X  ,j@»X  ,T)|dT  . 

1  (j-2)0  X  3  J*1  3  3“1 


In  order  to  estimate  the  last  integral  recall  that 


k 

% 


|Nx(xj<P;xj_1,T) I  <  Kg(P  -  T)"1/2,  p  >  T  . 

Substituting  this  in  the  estimate  of  J3  we  obtain  after  some  algebra 

U3I  <  k/3  . 

Analogous  considerations  yield  |J£|  ♦  |J4I  +  |JJ|  <  Kg/0.  Therefore  the  error  terns 
e^tt)  can  be  estimated  by 

l*g(t) |  <  Kg/I  +  K10(j  -  1)'3/20_V2  . 

For  the  error  Eg  in  <4. 6)  we  have 

p-t  j9  t 

|E»|  <  l  I  |e«(T> |dT  +  /  |e|(T)|dx  < 

j-1  (j-1)0  (p-1)« 

9  .  p  j9 

<  /  |eI(T)|dx  +  l  /  |ei(T)|dr  . 

0  8  j-2  (J-1)6 


Now  trivially  |e^(T)|  <  2C? . 


Therefore 


|Bel  <  2C  6  +  K  /5  f  8  +  K  /5  !  (j  -  1>‘3/2  . 
j-2  j-2 

Since  f  9  <  p6  <  T,  and  the  series  I  (j-1)  3/>2  converges,  we  obtain 
j-2  j>2 

(4.11)  |Eel  <  K^/5  . 

This  proves  the  lemma  • 
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5.  The  error  estimate i 

The  purpose  of  this  section  is  to  give  an  estimate  of  the  speed  of  convergence  of 
the  approximate  interface  Sg(t)  to  the  true  interface  s(t)  .  In  view  of  Lemma  4.1 
and  (4.1)  this  will  also  complete  the  proof  of  the  theorem. 

Lemma  5,1 1  There  exists  a  constant  C  depending  upon  H,  b,  Gj,  G2/  T  such  that 

,a9  '  a,T  *  C/S  * 

Proof:  Subtract  (4.6)  from  (4.2)  and  use  (4.11)  to  obtain 


|.(t)  -  *e(t)|  <  /  |g(u(0,T) ,T)  -  g(ue(0,TJ,T)|<Jt  + 


s(t)  s9(t) 

+  |/  u(x,t)dx  -  /  u8(x,t)dx|  + 

0  0 

sft(t)  X 

o  P 

+  |  /  uQ(x,t)dx  -  /  Ug(x,t)  dx|  +  , 


where  Ug  is  the  unique  solution  of  the  problem 

LUq  *0  in  P-  i  (0  <  X  <  8 g  ( t ) }  *  (0,T]  , 

6 

Ugx(0,t)  -  g(Ug(0,t) ,t) ,  0  <  t  <  T  , 

Ug(x,0)  -  h(x),  x  e  (0  ,b]  , 

Ug (8g (t) / 1 )  -  0,  0  <  t  <  T  . 

On  the  basis  of  the  maximum  principle  and  standard  barriers  estimates  £10,71,  as  in 

Lemma  3.3,  we  have 


0  <  u(x#t)  <  C^stt)  -  x) , 


(x,t)  e  dt  , 


0  <  Ug(X,t)  <  ( s  g  ( t )  -  x) t 


(x,t)  e  P-  . 

*8 


Proposition  5.1:  There  exist  constants  and  B2  depending  upon  H,  b ,  G1#  G2, 

T ,  such  that 
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«(t>  Vfc>  _  t  "*  ,8(P)  -  “9(p)l 

j /  u(x,t)dx  -  /  u6(x,t)<ix|  <  B  /  - - - - 

0  0  0  /t  -t 

t  »«IT 

+  B,  /  -  —  dr  . 


0  -  t 


Proof  of  Proposition  5.1?  Set 

a(t)  «  min{s<t) ,s9(t))>  B(t)  -  max(s(t) ,s0(t) } , 

X(t)  -  5(t)  -  5(t) ,  t  e  [0,t]  . 

Obviously  5(M,  5(*)  are  non- decreasing  Lipschitz  continuous  functions  with 
Lipschitz  constant  bounded  by  C-.  Then 


s ( t )  9  m  oi(t) 

|  /  U(  X  r  t )  dX  -  /  Ue(xft)dx|  <  J  |u(Xj,t)  -  ufi(x,t)|dx  + 

0  0  0 


3(t) 

/  |y(x,t) | dx  *  +  i2 

«(t) 


y(x,t) 


u(x,t)  if  o(t)  “  8^(0 

u^(x,t)  if  5(t)  -  s(t) 


He  dominate  the  integrand  in  1^  by  the  sum  |v  -  u^l  +  |u  -  v|  where  v  is  defined 
in  Lemma  3*3  and  v^f  v2  are  the  solutions  of  the  problems 

[  Lv«  -  0,  {0  <  x  <  S(t)}  *  (0  <  t  <  T>  , 


Lvt  -  0, 

vix(0.t)  -  0, 

Vj (x,0)  *  0, 


0  <  t  <  T  , 
0  <  x  <  b  , 


Vj(a(t)#t)  «  0  <  t  <  T 


Lv2  -  0,  (0,*)  x  (o,Tj  , 

v2x<0,t)  -  -|g<u(0,t),t)  -  g<u9(0,t)ft)|f  0  <  t  <  T  , 

t  x  10  • 


v2  ( x ,  0 )  *  0 


As  for  v2#  it  can  be  represented  explicitly  [9],  by 

v2(x,t)  -/  N(x,ti0#T)|g(u(0,T),T)  -  g(ufl(0,T) ,T) |dt  . 

0 

Therefore  by  virtue  of  the  Lipschitz  continuity  of  g(*,t)  and  (4.1)  we  obtain 

G1  ,*  ,4,x 

v2(xrt)  <  —  C2  /  ■  ■■■  ■—  dt  . 

/*  *  0  /t  -  i 

By  an  argument  of  page  87,  v^  can  be  dominated  by  z(x,t)  +  z(-x,t)  where 
solves 

Lz  «  0,  (-a(t)  <  x  <  •}  x  {0  <  t  <  T>  , 

(Pz)  •  *(-a( t),t)  -  C,5(t)  ,  0  <  t  <  T  , 

z(x,0)  *  0,  -b  <  x  <  m  . 

Then  there  exists  a  constant  K12  depending  only  upon  the  data  such  that 

t 

z{x,t)  <  K  /  iSl  |T  (x#t|-u(T) ,x) | dt  . 


Therefore 


o(t)  a(t)  •  t 

/  v  (x,t)dx  <  /  z (x,t) dx  <  /  z(x,t) dx  <  K  / 

0  1  -  -  13  0 
-o(t)  -a(t) 


dt  . 


The  estimate  of  I2  is  done  by  exploiting  again  the  methods  of  (2)  page  07  and 
dominating  the  integrand  by  z(x,t),  the  solution  of 

Lz  «  0,  {5<t>  <  x  <  *}  *  (O  <  t  <  t)  , 

(P-)  z(x,0)  -  0,  b  <  x  <  •  , 

k  z(a(t),t)  -  C1?(t),  0  <  t  <  T  . 


It  gives 


t  Ijl 

*2  ‘  *14  /  7=  dT  * 

0  /t  -  T 


Proposition  5.2 t  There  exist  a  constant  B  depending  upon  H,  b,  G^,  G^,  T  such 
that 


8e(t)  XP 

l/  Ug(x,t)dx  -  /  U0(x,t)dx|  <  B0  . 


Proof  of  Proposition  5.2:  For  t  >  0  fixed  we  have  s0(t)  >  •  Then 


s0(t)  x  x 

9  p  p 

I /  uQ(x,t)dx  -  /  uQ<x,t)dxj  <  /  |uQ(x,t)  -  u0(x,t)|dx 

0  0  0 

S0(t) 


+  /  uA(x,t)dx  -  ♦  Jn  . 


1  2 


Standard  barrier  arguments  [7,8,9]  give 

sup  | u Q ( x , t ) |  <  C^(b  +  C^T)  . 

(x,t)e£>- 

s0 

On  the  other  hand  by  the  construction  of  sfl(*),  the  distance  sa(t)  -  x  is  at  most 

O  O  p 

Hence 

J2  4  *159  * 

To  estimate  the  integrand  in  J-j  we  proceed  as  in  Lemma  4*1.  In  fact,  in  the  present 

situation  the  estimates  are  simpler  since  u0  -  u0  has  zero  flux  at  the  fixed  face 

x  -  0.  We  obtain  the  estimate  analogous  to  < 4 • 1 ) 

|uA(x,t)  -  U a( X  f t)  |  <  K  sup  I S A ( T )  -  Sfl(T)|  . 

0  0  16  0<x<t  9  9 


As  remarked  above  for  all  t  e  [0,T]  we  have 

se(t)  *  se(t)  <  C^B  9 
and  hence  the  proposition  is  proved. 
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We  now  conclude  the  proof  of  Lemma  5.1.  Notice  that  virtue  of  the  Llpschltz 


continuity  of  g(*,t)  and  (4.1)  we  have 


/  lg(u(0,T) ,T)  -  g(u.(0 

0 


t,  t  |u(0,T)  -  Ufl(0,T) I 
/  T  )  ,  T  )  |  dT  -  '  - 


di 


t  -  T 


V  *  I<5It 
<  G  T  /2C  /  ■  ■  ■■  <JT  . 

0  /t  -  T 


Putting  together  the  various  estimates  so  obtained  we  see  that  for  all  0  <  t  <  T 
have 


t  iSi  t  i«i 

\a(t)  -  ;e(t)  |  <  K  /e  +  K  /  -  ■  -  dT  +  K  7  /  - - dT 

0  mT-T  0  /t  -  T 


Now 


|s(t)  -  80(t)|  >  |s(  t)  -  SQ(t)|  -  1 S  0(  t )  -  se(t)| 


and 


l?l  -  sup  |s(p)  -  8a(P)|  <  ISI  +  sup  |sfl(P)  -  8a(P) |  <  HI  +  C,6  . 
0<P<T  T  0<P<T  8  8  T  1 


Consequently  the  above  implies 


»«»  <  K,„/e 

t  20 


+  K. 


21 


/ 


0 


dT  . 


The  proof  of  the  lemma  is  concluded  with  an  application  of  Gronwall's  inequality. 


we 
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6.  Modifications  of  the  achgme: 


Consider  the  sequence  of  approximating  problem  (P^)  introduced  in  Section  2.  The 
theorem  remains  true  if  we  modify  the  flux  condition  on  the  fixed  face  x  *  0  in  (P^) 
as  follows 

(6.i)  Retarded  flux. 

We  set  u£(0,t)  »  g<h(0),t),  0  <  t  <  8  and  for  j  >  1  we  replace  the  flux 
condition  in  (Pj)  with 

u£<0,t>  *  g(u^"**(0,t  -  8),t),  (j  -  1)8  <  t  <  j8  . 

(6.ii)  Piecewise  constant  flux. 

The  modification  oonsists  in  freezing  the  flux  on  x  •  0  in  the  rectangle 

at  its  value  at  the  lower  left  corner . 

Set  u^(0,t)  «  g(h(0),0),  0  <  t  <  8  and  then  for  j  >  1 

ul(0,t)  -  -  1)6), (j  -  1)9),  (j  -  1)0  <  t  <  j8  . 

The  proof  of  the  convergence  and  the  error  estimate  for  the  scheme  obtained  with  the 
above  modifications,  is  carried  out  in  essentially  the  same  way  as  indicated  in  the 
previous  sections.  Some  minor  modifications  are  needed  which  we  leave  to  the  reader, 
(cf.  [1]). 


j 

1 

i 
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